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We investigate the properties of the FLRW fiat cosmological models in which the vacuum energy 
density evolves with time, A(t). Using different versions of the A(t) model, namely quantum field 
vacuum, power series vacuum and power law vacuum, we find that the main cosmological functions 
such as the scale factor of the universe, the Hubble expansion rate H and the energy densities are 
defined analytically. Performing a joint likelihood analysis of the recent supernovae type la data, 
the Cosmic Microwave Background (CMB) shift parameter and the Baryonic Acoustic Oscillations 
(BAOs) traced by the Sloan Digital Sky Survey (SDSS) galaxies, we put tight constraints on the 
main cosmological parameters of the A(t) scenarios. Furthermore, we study the linear matter 
fluctuation field of the above vacuum models. We find that the patterns of the power series vacuum 
A = m H + n2 H 2 predict stronger small scale dynamics, which implies a faster growth rate of 
perturbations with respect to the other two vacuum cases (quantum field and power law), despite 
the fact that all the cosmological models share the same equation of state (EOS) parameter. In 
the case of the quantum field vacuum A = no + n,2 H 2 , the corresponding matter fluctuation field 
resembles that of the traditional A cosmology. The power law vacuum (A oc a~ n ) mimics the 
classical quintessence cosmology, the best fit being tilted in the phantom phase. In this framework, 
we compare the observed growth rate of clustering measured from the optical galaxies with those 
predicted by the current A(t) models. Performing a Kolmogorov-Smirnov (KS) statistical test 
we show that the cosmological models which contain a constant vacuum (ACDM), quantum field 
vacuum and power law vacuum provide growth rates that match well with the observed growth 
rate. However, this is not the case for the power series vacuum models (in particular, the frequently 
adduced A oc H model) in which clusters form at significantly earlier times (z > 4) with respect to all 
other models (z ~ 2). Finally, we derived the theoretically predicted dark-matter halo mass function 
and the corresponding distribution of cluster-size halos for all the models studied. Their expected 
redshift distribution indicates that it will be difficult to distinguish the closely resembling models 
(constant vacuum, quantum field and power-law vacuum), using realistic future X-ray surveys of 
cluster abundances. However, cluster surveys based on the Sunayev-Zeldovich detection method 
give some hope to distinguish the closely resembling models at high redshifts. 
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1. INTRODUCTION 

Over the past decade, studies of the available high 
quality cosmological data (supernovae type la, CMB, 
galaxy clustering, etc.) have converged towards a cos- 
mic expansion history that involves a spatially flat ge- 
ometry and a recent accelerating expansion of the uni- 
verse (cf. [H, [2, !, i, i, H and references therein). 
This expansion has been attributed to an energy com- 
ponent called dark energy (DE) with negative pressure, 
which dominates the universe at late times. The sim- 
plest type of DE corresponds to the cosmological con- 
stant (see for reviews). The so-called concordance 
model (or ACDM model), which contains cold dark mat- 
ter (DM) to explain clustering, flat spatial geometry and 
a cosmological constant, A, fits accurately the current ob- 
servational data and thus it is an excellent candidate to 
be the model that describes the observed universe. How- 
ever, the concordance model suffers from, among others 
two fundamental problems: (a) The "old" cosmo- 



logical constant problem (or fine tuning problem) i.e., the 
fact that the observed value of the vacuum energy den- 
sity (pa = Ac 2 /87rG ~ 10~ 47 GeV 4 ) is many orders of 
magnitude below the value found using quantum field 
theory (QFT) 0], and (b) the coincidence vroblem [IlT] 
i.e., the fact that the matter energy density and the vac- 
uum energy density are of the same order just prior to 
the present epoch, despite the fact that the former is a 
rapidly decreasing function of time while the latter is just 
stationary. The extremal possibility concerning prob- 
lem (a) occurs when the Planck mass Mp ~ 10 19 GeV 
is used as the fundamental scale; then the ratio Mp/p\ 
becomes ~ 10 123 . One may think that physics at the 
Planck scale is not well under control and that this enor- 
mous ratio might be fictitious. However, let us consider 
the more modest scale v = 2Mw/g — 250 GeV associ- 
ated to the electroweak Standard Model (SM) of Particle 
Physics (the experimentally most successful QFT known 
to date), where Mw and g are the boson mass and 
the SU(2) gauge coupling, respectively. In this case, the 



2 



ratio of the predicted vacuum energy versus the mea- 
sured one is still very large: it reads \(V)\/pa > 10 55 , 
where (V) = —(l/8)MjjV 2 is the vacuum expectation 
value of the Higgs potential and Mjj > 114.4 GeV is the 
lower bound on the Higgs boson mass in the SM. 

Attempts to solve the above problems have been pre- 
sented in the literature (see [1, [t| [l2j and references 
therein) . Historically, the first attempts to solve the "old 
A problem" (the toughest one in the list) were based on 
dynamical adjustment mechanisms [l3| aiming to avoid a 
mere fine tuning of the various QFT contributions. Later 
on it was realized that a dynamical A could also be useful 
to overpass the coincidence problem. The idea is to re- 
place the constant vacuum energy either with a DE that 
evolves with time (quintessence and the likeQ), or al- 
ternativel y w ith a time varying vacuum energy density, 
PA(t) (cf. [14l. [l5l|). In the quintessence context, one can 
introduce an ad hoc tracker scalar field cj}^ (different 
from the usual SM Higgs field) rolling down the poten- 
tial energy V{4>), and therefore it could resemble the DE. 
Detailed analysis of these models exist in the literature, 
including their confrontation with the data fl7l. [i~8l. fl9L [20| 
- see also [1, H, [2l[ for extensive reviews. Nevertheless, 
the idea of a scalar field rolling down some suitable poten- 
tial does not really solve the problem because the initial 
value of the DE still needs to be fine tuned. Moreover, 
the typical mass of <f> is usually very small. Indeed, being 
4> unrelated to the SM physics, it is usually thought of 
as a high energy field (0) ~ Mx where Mx is some scale 
near the Planck mass. If one assumes the simplest form 
for its potential, namely V(<p) — m|0 2 /2, and requires 
it to reproduce the present value of the vacuum energy 
density, p\ — (V((f>)) ~ 10~ n eV 4 , the mass of <j> is pre- 
dicted in the ballpark of ~ H ~ 10~ 33 eV. This 
is an inconceivably small mass scale in Particle Physics. 
Thus, the problem that one is creating along with the 
introduction of <fr is much bigger than the problem that 
one is intending to solve, for one is postulating a mass 
scale which is 30 orders of magnitude below the mass 
scale 771a associated to the value of the vacuum energy 
density m\ = p 1 / 4 ~ 2.3 x 10~ 3 eV. 

Current observations do not rule out the possibility of 
a dynamical DE 0, @, i, 0, H, @]. They indicate that 
the dark energy EOS parameter w = Pde/pde is close 
to —1 to within ±10%, if it is assumed to be constant, 
whilst it is much more poorly constrained if it varies with 
time [2j . Let us note, interestingly enough, that models 
with running vacuum energy, may appear in practice with 
a non-trivial "effective EOS" to — w(t) [13], which can be 
accommodated by observations. Indeed, the basic fea- 
ture in this cosmological ideology is that, although the 
EOS parameter of the vacuum energy is strictly equal 
to -1, a time evolving A generally results in an effective 
EOS, usually a function w = w(a) of the scale factor that 
approaches to — > —1 at the present time. This effective 
EOS may be the result either of the fact that we are 
trying to describe the variable A(i) model as if it were a 
quintessence type model with conserved DE [22j , or of the 



existence of other DE entities mixed up with the variable 
A(i), as e.g. in the case of the AXCDM model [H HI- 

Although the precise functional form of A(i) is still 
missing, an interesting QFT approach within the con- 
text of the renormalization group (RG) was proposed 
long ago [24], H^]. Later on, the RG-running framework 
was further explored in [2(| from the viewpoint of QFT in 
curved space-time by employing the standard perturba- 
tive RG-techniques of Particle Physics (see also [27j for a 
functional RG approach in a non-perturbative context). 
A more phenomenological point of view was addressed in 
[28| . in which a time varying A(t) parameter could be a 
possible candidate to solve the two fundamental cosmo- 
logical puzzles. There is an extensive literature on time 
varying A models (cf. [29|, [3(1 HH and references therein). 

It is worth noticing, that for an important class of A(t) 
models there is a coupling between the time-dependent 
vacuum and matter. A first work confronting supernovae 
data with a (RG- inspired 26]) model of this kind was pre- 
sented in [32| - see also J33| . Using the combination of 
the conservation of the total energy with the variation of 
the vacuum energy, one can show that these A(t) models 
provide either a particle production process or that the 
mass of the dark matter particles increases. Despite the 
fact that most of the recent papers on these matters are 
based on the assumption that the DE evolves indepen- 
dently of the dark matter, the unknown nature of the DE 
and DM implies that at the moment we can not exclude 
the possibility of interactions in the dark sector, whether 
at the level of a variable A(t) [H, HH or from coupled 
quintessence [lij . Another possibility is that matter is 
strictly conserved and that both A(t) and the gravita- 
tional coupling G(t) are running [35|, [36[, but this option 
will not be scrutinized here. 

The scope of the present work is to study the obser- 
vational consequences of the overall dynamics for a wide 
class of time varying vacuum energy models, in interac- 
tion with matter, in the light of the most recent cosmo- 
logical data. The structure of the paper is as follows. 
The basic theoretical elements of the problem are pre- 
sented in section 2, where we introduce [for a spatially 
flat Friedmann-Lemaitre-Robertson- Walker (FLRW) ge- 
ometry] the basic cosmological equations. In sections 3 
and 4 we place constraints on the main parameters of our 
vacuum models by performing a joint likelihood analysis 
utilizing the so called 'Constitution set' of the SNIa data 
@, the shift parameter of the Cosmic Microwave Back- 
ground (CMB) and the observed Baryonic Acoustic 
Oscillations (BAOs; [13] )• Also, we investigate whether 
the current vacuum cosmological models can yield a late 
accelerated phase of the cosmic expansion. In section 5 
we compare the observed linear growth rate of clustering 
measured from the SDSS galaxies with those predicted 
by the A(t) models, explored here, while in section 6 we 
present the corresponding theoretical predictions regard- 
ing the formation of the galaxy clusters and the evolu- 
tion of their abundance. In section 7 we draw our conclu- 
sions. Note, that throughout the paper we use Hq = 70.5 
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km/sec/Mpc [EH^. 



2. COSMOLOGY WITH A TIME DEPENDENT 
VACUUM 

In the framework of a spatially flat FLRW geometry 
the basic equations which govern the global dynamics of 
the universe are 



Pa = 3H 2 



(1) 



and 



p'm + PA + 3H(p m + P m + p A + P A ) = , (2) 



where p m is the matter energy density and pa is the 
vacuum energy density, while P m and Pa are the corre- 
sponding matter and vacuum pressures. Notice, that in 
order to simplify our formalism we use geometrical units 
(8ttG = c = 1) in which p A = A. 

In the present work, we would like to investigate the 
potential of a time varying A = A(t) parameter to ac- 
count for the observed acceleration of the expansion of 
the Universe. Within this framework it is interesting to 
mention that the equation of state takes the usual form 
Ph(t) = —pA(t) = — A(t), and we remark that this EOS 
does not depend on whether A is strictly constant or 
variable. Also, by introducing in the global dynamics 
the idea of the time-dependent vacuum, it is possible to 
explain the physical properties of the DE as well as to 
ameliorate the status of the fine tuning and the coin- 
cidence problem respectively. In the matter dominated 
epoch (P m = 0), eq.© leads to the following useful for- 
mula: 



Pn 



3H Pm = - A, 



and considering eq. (pj, we find: 



H 



-H z 



A 

2' 



(3) 



(4) 



where the over-dot denotes derivatives with respect to the 
cosmic time. If the vacuum term is negligible, A{t) — > 0, 
then the solution of the above equation reduces to that of 
the Einstein de-Sitter model, H(t) = 2/3i, as it should. 
Obviously, if we consider the case of A(t) ^ 0, then it 
becomes evident - see eq.® - that there is a coupling 
between the time-dependent vacuum and matter compo- 
nent. This equation was first considered by M. Bron- 
stein in a rather old paper [39(. Note, that the tradi- 
tional A = const cosmology (or ACDM model) can be 
described directly by the integration of the eq.© [for 
more details see section 4.1], but this same equation is 
also valid for A = A(t), in which case a supplementary 
equation for the time evolution of A is needed in order 
to unveil the dynamics of this model. As we have al- 
ready stated before, the link in eq.© between p' m and A 



is very important because interactions between DM and 
DE could provide possible solutions to the cosmological 
coincidence problem. This is the reason for which sev- 
eral papers have been published recently in this area (34| 
proposing that the DE and DM could be coupled. Alter- 
natively, one may use a variable A together with another 
entity X, that ensures a total self-conserved DE density, 
Pde = PA + Px, such that the ratio pde/ Pm remains 
fairly stable for most of the Universe's lifetime fl4j . In 
both cases a time-dependent DE is needed. 

In the following, we are going to concentrate on models 
where the time dependence of A appears always at the 
expense of an interaction with matter. In this context, 
the corresponding time evolution equation for the matter 
density contrast D = 5p m / p m , in a pressureless fluid, is 
given by : 



D + (2H + Q)D 
where 



- 2HQ - Q 



D = Q. 



3H 2 -A Q(t) = -A/p„ 



(5) 



(6) 



It becomes clear, that the interacting vacuum energy af- 
fects the growth factor via the function Q{t). Of course, 
in order to solve the above differential equations we need 
to define explicitly the functional form of the A(t) com- 
ponent. Notice that the approach based on eq. © effec- 
tively implies that the DE perturbations are negligible 
(i.e., in this case we set SA = 0). This is justified in most 
cases, specially for perturbations well inside the sound 
horizon (of the dark energy medium) where it behaves 
very smoothly [4lL l42l] . In fact, one can explicitly derive 
eq. ([5]) starting from the general coupled set of matter and 
dark energy perturbations [see equations (17), (25) and 
(27) of Ref. (4l|, which we refrain from repeating here]. 
Assuming that matter and A(t) interact as in eq.©, tak- 
ing the limit where the DE perturbations are neglected, 
and assuming also that the produced particles of matter 
have negligible velocities with respect to the co-moving 
observers, then the aforesaid coupled set of matter and 
DE perturbation equations lead to a second order differ- 
ential equation for D that boils down to eq. ([5]) above. 
In this way, we can state that this effective equation fol- 
lows from the general relativistic treatment of perturba- 
tions, within the aforementioned set of approximations, 
whereas in [401 ] it was originally proven directly within 
the Newtonian formalism. 



3. LIKELIHOOD ANALYSIS 

In this work, we use a variety of cosmologically relevant 
observations in order to constrain the vacuum models 
explored here (see section 4). These are: 

• Baryonic Acoustic Oscillations (BAOs): 

These are produced by pressure (acoustic) waves 
in the photon-baryon plasma in the early universe, 
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generated by dark matter overdensities. Evidence 
of this excess was recently found in the clustering 
properties of the luminous SDSS red-galaxies [37| 
and it can provide a "standard ruler" with which 
we can constraint the dark energy models. In par- 
ticular, we use the following estimator 



A(p) = 



z 2 E(a s )]^ 



da 



a 2 E(a) 



2/3 



(7) 



measured from the SDSS data to be A = 0.469 ± 
0.017, with z s = 0.35 [or a s = (1 + Zs)- 1 ~ 0.75] 
and E(a) = H(cl)/Hq is the normalized Hubble 
flow. Therefore, the corresponding Xbao function 
is simply written 



Xbao(p) 



[A(p) - 0.469] : 
0.017 2 



(8) 



where p is a vector containing the cosmological pa- 
rameters that we want to fit. 

• CMB shift parameter: A very accurate and 
deep geometrical probe of dark energy is the angu- 
lar scale of the sound horizon at the last scattering 
surface as encoded in the location if T of the first 
peak of the Cosmic Microwave Background (CMB) 
temperature perturbation spectrum. This probe, 
called CMB shift parameter 0, Q EH , is defined 
as: 



R 




da 



a 2 E{a) 



(9) 



Note, that the Hubble function H(a) includes also 
the radiation component (fl r ~ 10~ 4 ). The shift 
arameter measured from the WMAP 5-years data 
is R = 1.71 ± 0.019 at z is = 1090 [or a Ls = 
(1 + zis)" 1 ~ 9.17 x 10~ 4 ]. In this case, the xLb 
function is given by 



xLb(p) = 



[R(p) - 1-71] 2 



(10) 



0.019 2 

Note that the measured CMB shift parameter is 
somewhat model dependent but mostly to models 
which are not included in our analysis. For ex- 
ample, in the case where massive neutrinos are in- 
cluded or when there is astrongly varying equation 
of state parameter. The robustness of the shift pa- 
rameter was tested and discussed in [46]. 

• SNIa distance modulii: We additionally utilize 
the 'Constitution set' of 397 type la supernovae 
of Hicken et al. [f|. In order to avoid possible 
problems related with the local bulk flow, we use a 
subsample of the overall sample in which we select 
those SNIa with z > 0.023. This subsample con- 
tains 351 entries. The corresponding XsNia function 
is: 



XSNIa(P) = 



351 

E 



/* tb (ai,p)-M ob "(ai) 



(11) 



where aj = (1 + Zi)^ 1 is the observed scale factor 
of the Universe, Zi is the observed redshift, /i is the 
distance modulus corresponding to flat space: 

H = m- M = 51ogrf L + 25 (12) 

and <iz,(a,p) is the luminosity distance 

d L (a,p) = - t-Mri, (13) 

aJ a y H(y) 

Note, that c is the speed of light (c = 1 here) . 

We can combine the above cosmologically tests, using a 
joint likelihood analysis, in order to put even more strin- 
gent constraints on the free-parameter space, according 
to: 

Aoi(p) = £bAO x £cmb X -CsNIa 



or 



xL(p) = Xbao 



Xcmb 



XSNIa 



(14) 



with the likelihood estimator defined as Cj oc 
exp[— Xj/2]. Note, that we sample the Q m parameter in 
steps of 0.01 in the range [0.1,1] and that the likelihoods 
are normalized to their maximum values. We will report 
la uncertainties of the fitted parameters. Note that the 
overall number of data points used is N to t = 353 and the 
degrees of freedom: dof— Ntot ~ ^fit; with net the number 
of fitted parameters, which vary for the different models. 



4. CONSTRAINTS ON THE TIME EVOLVING 
VACUUM MODELS IN FLAT SPACE 

As we have already mentioned, the exact nature of a 
possible time varying vacuum has yet to be found. A 
large number of different phenomenological parameteri- 
zations have appeared in the literature treating the time- 
dependent A(i) function. For example, the authors of [29j 
considered A(t) = 3jH 2 , with the constant 7 being the 
ratio of the vacuum to the sum of vacuum and matter 
density (see also [13]), while [28[ and [13] proposed a 
different ansatz in which A(t) oc a~ 2 . Also, several pa- 
pers, (see for example [H, [48| and references therein) 
have investigated the global dynamical properties of the 
universe considering that the vacuum energy density de- 
creases linearly with the matter energy density, A cx p m . 
Carneiro et al. [13, HFJ] used a different pattern in which 
the vacuum term is proportional to the Hubble param- 
eter, A(a) oc H(a), while [HH considered a power series 
form in H . Attempts to provide a theoretical explana- 
tion for a dynamical A(t) have also been presented in 
the literature using the renormalization group (RG) in 
quantum field theory (see [TJ, [13, US H3] and references 
therein). The RG- inspired form for a QFT running vac- 
uum in curved space-time is 
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A(H) = n + n 2 H 2 



(15) 



where no and n 2 are constants. This evolution law can 
mimic the quintessence or phantom behavior |53| as well 
as a smoothly transition between the two [221]. Notice, 
that the functional form of eq. (fT5|) has been also used in 
[l5| . Such a form is indeed crucially different from just 
considering that the vacuum energy is proportional to 
H 2 , in the sense that the former is an "affine quadratic 
law" (n 7^ 0). Remarkably, the structure of eq. p3|) can 
be motivated from the QFT framework of anomalous in- 
duced inflation (cf. [36| and references therein). In an- 
other vein, the aforementioned possibility that the vac- 
uum energy could be evolving linearly with H has been 
motivated theoretically in the literature through a possi- 
ble connection of cosmology with the QCD scale of strong 
interactions [53j . This option, however, is not what one 
would expect from re-normalizable QFT in curved space- 
time because from general covariance we should rather 
expect even powers of H, as e.g. in the law (fTS)). There 
is, however, the possibility to add non-analytic terms in 
the effective action (see [54( for a recent and interesting 
attempt in this direction). Such a linear dependence in 
H has also been proposed from a possible link of DE with 
QCD and the topological structure of the universe [5^ . 
Let us, however, note that a connection with QCD can 
also be achieved through a relaxation mechanism of A 
and without using the hypothesis of linearity in H - see 

In this work, we consider a large family of flat vacuum 
models and with the aid of the current observational data 
(see section 2) we attempt to put stringent constraints on 
their free parameters. Also, we investigate thoroughly 
the time evolution equation of the mass density contrast 
in the linear regime as well as the formation of galaxy 
clusters and the evolution of their abundance. In the fol- 
lowing subsections, we briefly present these cosmological 
models which trace differently the vacuum component. 



4.1. The standard A-Cosmology 

Without wanting to appear too pedagogical, we remind 
the reader of some basic elements of the concordance A- 
cosmology in order to appreciate the differences with the 
A(t) models explored subsequently. The vacuum term in 
eq.(d]) is constant and given by A = 3HaHo- Therefore, 
it is a routine to integrate eq.((4]) and obtain the Hubble 
function: 



H(t) 



Ho coth 



(16) 



where fl\ = 1 — f2 m . Note, that £l m is the matter density 
parameter at the present time. Using now the definition 
of the Hubble parameter H = a /a, the scale factor of the 



universe a(t), normalized to unity at the present epoch, 
evolves with time as: 



a(t) 



£7 m 



1/3 



sinhl (**°^ t 



The cosmic time is related with the scale factor as 



t(a) 



-sinh 




(17) 



(18) 



Combining the above equations we can define the nor- 
malized Hubble expansion as a function of the scale fac- 
tor: 



E 2 (a) 



H 2 (a) 



= VL\ + f2 m a 



-3 



(19) 



The inflection point [namely, the point where the Hubble 
expansion changes from the decelerating to the acceler- 
ating regime, d(ti) = 0] takes place at: 



ti 



-sinh 








j , a>i = 


2n A 



"I 1/3 



(20) 



Comparing the concordance model with the observa- 
tional data we find that the best fit value, within the lex 
uncertainty, is O m = 0.28 ± 0.01 with xL(^m) ^ 431.2 
(<io/=352Y, which is in good agreement with recent stud- 
ies [l|, H, H, H, @] . Therefore, the current age of the uni- 
verse is to — 13.9Gyr while the inflection point is located 
at t] ~ 0.52io, a] ~ 0.58 (hence at redshift Z] ~ 0.72). 
Let us mention that some recent (approximately model- 
independent) determinations of this point suggest it to lie 
at a more recent time (lower redshift) [56| . although the 
results are still compatible with the ACDM value within 
2cr. 

Finally, solving eq.© for the A cosmology [Q(t) = 
—A/p m = 0], we derive the well known perturbation 
growth factor (see [57|): 



D{a) 



dx 



o x 3 E 3 (x) 



(21) 



In particular, for E(a) = flm 2 a~ 3 / 2 it gives the stan- 
dard result D(a) = a, which corresponds to the matter 
dominated epoch, as expected. Notice, that (|2"Tj) is nor- 
malized to unity at the present time, to, because in our 
convention a(to) = 1 - cf. equations (| 1 71 1 8[) . 



4.2. The A(t) model from quantum field theory 

Let us consider the vacuum solution (fT5|) proposed 
in [2l,[3l,[3l,[3l| using the renormalization group (RG) in 
quantum field theory (hereafter Ana model). We select 
the coefficients in that equation as no = 3 Hq (J1a — 7) 
and ri2 = 37. Therefore, 
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A(H) = A + 3 7 (H 2 - Hi) . 



(22) 



In this way, the vacuum energy density is normalized 
to the present value: A = A(H ) = 3Qa#o- Without 
going into the details of that model, let us recall that 
7 (called v in the above papers) is interpreted in the 
RG framework as a "/3-function" of QFT in curved-space 
time, which determines the running of the cosmological 
constant. The predicted value is 
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a M 2 
12tT Mj 



(23) 



where Mp is the Planck mass and M is an effective 
mass parameter representing the average mass of the 
heavy particles of the Grand Unified Theory (GUT) near 
the Planck scale, after taking into account their multi- 
plicities. Since a = ±1 (depending on whether bosons 
or fermions dominate in the loop contributions), the co- 
efficient 7 can be positive or negative, but | 7 | is natu- 
rally predicted to be much less than one. For instance, 
if GUT fields with masses Mi near Mp do contribute, 
then |7| < 1/(12tt) ~ 2.6 x 1CT 2 , but we expect it to 
be lesser in practice because the usual GUT scales are 
not that close to Mp. By counting particle multiplici- 
ties in a typical GUT,a natural estimate is the range 
7 = 1CT 5 - 10" 3 (see [3^ for details). 

We will assume here that A evolves in the form (|2"2")l 
while it interacts with matter as in eq.©. Alternatively, 
one may use a variable A of the form (|22[) that interacts 
with a variable gravitational constant, G = G(t), such 
that matter is conserved (see [H, [36j]), but we shall not 
deal with this option here because we wish to consider 
only the class of variable A models in interaction with 
matter. 

It is important to emphasize that the main motiva- 
tion for the evolution law eq. (|2"2"|) stems from the general 
covariance of the effective action in QFT in curved space- 
time [H, [ill, Hil - for a review, see e.g. [58[. One expects 
that deviations from a strictly constant vacuum energy 
appear as a result of having a non-trivial external metric 
that describes an expanding FLRW background. Since 
the expansion rate of this background is H, we expect a 
power series in H. However, for a re-normalizable for- 
mulation of the effective action of the vacuum, only even 
powers of the expansion rate H can appear, the lead- 
ing correction being of 0(H 2 ). The next-to-leading term 
would be of 0(H 4 ), the subsequent one of 0{H 6 /M 2 ) 
etc. At the present time, all of the higher order correc- 
tions arc phcnomcnologically irrelevant compared to the 
first curvature correction ~ 0(H 2 ) in eq. ([22|) . Therefore, 
it is natural to take just the leading form, as in eq. (|22[) . 
See, however, [5^| for a possible effect of H 4 terms eval- 
uated at the electroweak crossover scale. 

It is interesting to point out that the Ape model can be 
used in different formulations where it helps to alleviate 



the cosmic coincidence problem [3, Il5j . It could also 
have a bearing on the fine tuning problem [l5j . 

From equations ^ and (f2"2")l we can easily derive the 
corresponding Hubble flow as a function of time: 



Hit) = H a J — — - coth 
V 1 ~7 



■iW(fiA-7)(l-7)* 



(24) 

The scale factor of the universe a(t), evolves with time 
as 



a(t) = a\ sinh 3 ^— 
where 

f 1 



-*W(fiA-7)(l-7)* 



/ Sl r , 



V^A - 7 



3(1-1) 



, (25) 



(26) 



Inverting eq. (|25p we determine the cosmic time as a func- 
tion of the scale factor: 



t(a) 



3iW(nA -7X1-7) 



sinh 1 



3(l-7)/2' 



(27) 

The age of the universe is simply given by this expression 
at a = 1, namely at the point to = t(a = I). The cos- 
mic time at the inflection point of the universe evolution 
[a(i/) = 0] is found to be 



ti 



3ffoV(JiA- 7)(l-7) 



sinh 



'1 - 37 



(28) 



and the corresponding value of the scale factor reads 



ai = 



(1 - 3 7 ) Sir, 



2(fi A - 7 ) 



(29) 



It becomes clear, that for 7 < 1/3 (this parameter is 
expected to be very small in QFT) and for the usual 
values of Sl/y (Sl\ > 7) the above inflection point exists. 

In practice the condition 7 < 1/3 is amply satisfied in 
the context of RG-inspired models in which the vacuum 
energy evolves as in eq . (I2"2"j) , the reason being that, in 
this QFT context, 7 appears from eq. (|23|) in which the 
highest scale M of the particle masses is expected to lie 
at (or below) the Planck scale. Therefore, from QFT we 
expect 7^1 and this is indeed what the comparison 
with experimental data confirms (see below). It is thus 
enlightening to expand (|2"9"|) in the limit 7 -C 1: 



1/3 



1-7 



3 V Sl m j2a 



(30) 



Since the coefficient 7 in this formula is positive (for 
the current values of the cosmological parameters), we 
see that the transition point from deceleration to accel- 
eration occurs earlier in time than in the concordance 



7 



ACDM model (|20|) . whereas for 7 < it occurs at a 
more recent time. 

Let us also consider the behavior of the matter and 
vacuum energy densities in this model as a function of the 
scale factor. Starting from the conservation law [see eq[3], 
and then trading the time derivatives for derivatives with 
respect to the scale factor, we obtain 



dp™. 
da 



3 dA 

-Pm = --T- 

a da 



(31) 



Using this equation in combination with §3§ and (|22p. 
we arrive at a simple differential equation for the matter 
density, 



dp m 3 

-j— + - (1 - 7) p m = , 
da a 

whose trivial integration yields 

Pm(«) = Pmoa~ 3(1 ~ 7) ■ 



(32) 



(33) 



Here p m o is the matter density at the present time (a 
1). Similarly, we find 



A(o) - A + 



7PmO 

1-7 



-3(1-7) 



(34) 



Substituting ea. (f25|) in the last two equations one may 
obtain the explicit time evolution of the matter and vac- 
uum energy densities, if desired. It is important to em- 
phasize from eq. (|33[) that the matter density does no 
longer evolve as p m (a) — p m oa~ 3 , as it presents a cor- 
rection in the exponent. This is due to the fact that 
matter is exchanging energy with the vacuum and this is 
reflected in the corresponding behavior of A (a) in ea. (|3"4")) . 
Substituting eq. (|34"]) in ea. (|22p . we immediately obtain 



E 2 {a) = 



f2 A - 7 



-3(1-7) 



= 1 + fi 



1 — 7 1 — 7 

-3(l-7)_l 



1 



'7 



(35) 



where in the second step we have used the cosmic sum 
rule for flat space f2 m + £1a = 1- Needless to say, eq. ([33)1 
can also be obtained by eliminating the cosmic time from 
equations ([24]) and ((25)) . as one can check. It is worth 
noting, that the normalized Hubble flow in this model 
[see eq. (|35[l ] can be viewed in a similar formulation as 
that of the concordance cosmology (for more details see 
the appendix A) as follows: 



E 2 {a) = tt t 



-3(1-7) 



where 



a 



1 



7 



(36) 



(37) 



As expected, for 7 — > (O m ~ O m ) all the above 
equations boil down to the canonical form within the 



concordance model - cf. previous section. Thus, the 
traditional A cosmology is a particular solution of the 
Arc model with 7 strictly equal to 0. 

Comparing the Aug model with the observational data 
(we sample 7 € [—1,0.3] in steps of 0.001) we find that 



the best fit value is fL 



0.28 



+0.02 
-0.01 



(or Q r , 



0.281) 

and 7 = 0.002 ± 0.001 with xf ot (O ro ,7) ^ 431.2 for 
dof= 351. We remark that the best fit value that we 
have obtained for 7 using the combined set of modern 
SNIa+BAO+CMB data becomes significantly smaller 
than the one obtained in the old analysis of Ref. [32| . 
where only a limited set of 54 supernovae data was em- 
ployed 1 . In fact, here we have been able to further re- 
strict the parameter 7 (called v in [32]) and push its 
value to its natural small range (viz. 7 ~ 10~ 3 or below) 
as expected from general QFT considerations [36| - see 
eq.p]). 

Being 7 small and positive, eq. (|30|) predicts that the 
transition point from deceleration to acceleration should 
be slightly earlier in time (hence at larger redshift) as 
compared to the standard ACDM case. 

It is interesting to point out that the small 7 value that 
we have obtained from the combined SNIa+BAO+CMB 
data is nicely compatible with the result obtained for 
this parameter from the analysis of the matter power 
spectrum of that model, performed in [6(| - see also [6l[ . 
Here we will analyze also the linear perturbation regime 
for the various models. However, rather than focusing 
on the power spectrum we will test the implications of 
this model on structure formation, through the study of 
the growth rate d In D(a)/d In a, the formation of galaxy 
clusters and the evolution of their abundances (see sec- 
tions 5 and 6). In all cases we need to find the linear 
matter fluctuation field D(a), which we shall determine 
under the assumption of vanishing A perturbations. As 
we have explained in section 2, this is justified for pertur- 
bations well inside the sound horizon and assuming that 
the produced matter particles have negligible velocities 
with respect to the co-moving observers. 

Therefore, we now proceed in an attempt to analyti- 
cally solve the differential eq.® in order to investigate 
the matter fluctuation field of the RG model (|22|) in the 
linear regime. To do so, we change variables from t to a 
new one according to the transformation 



coth 



Also from equations 
useful relations: 



■ iW(^A - 7) (1 - 7) * 



(38) 



and (|35|) we get the following 







tt A - 7 



E(a) 



v 2 -i = 



S 1 rt 



f2 A - 7 



-3/3 



(39) 



1 However, using the current SNIa data alone our best fit values 
are in agreement with those found by l32t . 
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where (3 = 1 — 7. Using (|38[) and (|39|) we find, after some 
algebra, that equation §5§ takes on the form 

3f3 2 (y 2 -l) 2 D" +f(P)y(y 2 -l)D' -2[g(P)y 2 -i;(f3)}D = 0, 

(40) 

where primes denote derivatives with respect to y, and 

/(/?) = 2/3(6/3 -5), 5 (/3) = (2-/3)(3/3-2), 
V(/3) = /3(4 - 3/3) . (41) 



In deriving eq. (|40|) we have substituted the various terms 
in ((5]) as a function of the new variable [see eq. ([39|) ] and 
the cosmological parameters. For instance, from equa- 
tions ©, (HH), ([32D and (J3HJ), we have 

p m = S^^a-t)^ 2 -!) 



= 3 7 #o 



9 



7 



1-7 



-7^(0 A - 7 )(1-^). 



(42) 



Factors of Hq drop at the end of the calculation. The dif- 
ferential equation (|4*0)l can be brought into the standard 
associated Legendre form by an appropriate transforma- 
tion. The growth factor solving this equation reads 



D(y)=C(y 2 -l) : 



-9/3 /I 13 1 3 
^ V + 2'2'3^ + 2'"^ 



1 



- 1 

where the quantity F is the hypergeometric function and 
C is a constant (for more details see the appendix B). 
Inserting eq. ([55|) into cq. (|4*3")) and using (|3T|) . we finally 
obtain the growth factor D{a) as a function of the scale 
parameter: 



9g-4 



1 13 1 



(44) 
(45) 



where 



Ci = c 



(4-9/9)/6/3 



As a consistency check we note that for 7 = (/3 = 1), the 
corresponding normalized to unity growth factor derived 
from eq. (l4~4"l) provides the same results as those of the 
concordance cosmology (see eq l2"Tj) . From this analysis, 
it becomes clear that the overall dynamics, predicted by 
the Arc model, extends nicely to that of the usual A 
cosmology and connects smoothly to it. For the analysis 
of other RG-inspircd phcnomcnological time varying A 
models, see Elf. 



4.3. The A(t) oc H 2 model 

We now consider that the vacuum energj 
cays as: A = 3jH 2 (hereafter A# x see [2i 



density de- 
40j). This 



model corresponds in setting n = in eq. (fTof . and 
therefore it can be derived as a particular case of the 
previous model. In this context, the basic cosmological 
equations become 



H(t) = 



3(1-7)* 



o(t) = 



3(1-7) 



3(1-7) 



and the normalized Hubble flow is 



E(a) 



-3(l- 7 )/2 



(46) 



(47) 



Note, that the constant 7 lies in the interval < 7 < 1, 
while the vacuum energy density remains constant ev- 
erywhere, which implies that Q m (a) = 1—7. This is 
consistent with the cosmic sum rule, since 7 = tt\ for 
this model. From eq. (|4"6"| it is obvious that this model 
has no inflection point. 

If we change the variables from t to a then the time 
evolution of the mass density contrast (see eqJS]) takes 
the following form 

a 2 D" + |o(l + 3 7 )L>' - ~(1 + 7)(1 - 3 7 )L> = (48) 
a general solution of which is 



D(a) = da 1 ^ + C 2 a-^ 1+ ^' 2 



(49) 



where Ci and C2 are the corresponding constants. No- 
tice, that a growing mode is present in this scenario if 
and only if 7 < 1/3, which implies that cosmic struc- 
tures cannot be formed via gravitational instability in 
Kh x models with 7 > 1/3. From a theoretical viewpoint, 
it is obvious that the Kh 1 pattern modifies the Einstein 
de-Sitter model. 

We find that the current model is unable to fit the com- 
bined observational data (SNIa+BAO+CMB). Indeed, 
we find that 7 = 0.64t°;° 2 with XsNia+BAO - 459 - 3 - In- 
terestingly, the addition of one more point (the CMB 
shift parameter) increases the overall likelihood func- 
tion by a factor of — 2.4 [xlotil) — 1104.8 with 7 = 
0.95 ± 0.01]. Using any of the fitted values of 7, it then 
follows from eq. (|46"l) that for this model the scale factor 
evolves as a ~ i e with e > 1. Under these conditions, 
this model is free from the horizon problem. However, 
unfortunately, it is unable to fit the present observational 
data, and moreover combining the latter statistical result 
with the 7 < 1/3 constrain we conclude that the current 
cosmological model is ruled out at high significance level. 



4.4. A power series A(t) model 

In this case, we parametrize the functional form of A(t) 
using a power series expansion in H up to the second 
order (see 51]) and assuming that there is no constant 
term: 



A = mH + n 2 H 2 



(50) 
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(hereafter Aps x ) . Performing the integration of eq. (|4|) we 
derive the following Hubble function: 



corresponds to a{tj) = and tj < to- Differentiating 
twice eq . ([5T>)) . we simply have: 



H(t) 



3 nit/2 



7 e' 



mt/2 



(51) 



where we have defined 7 = 3 — n2 and expect 7 to remain 
around 3 (or |na| <C 1). In fact, we do not foresee that 
the coefficient could be large (for similar reasons as in 
section 4.2). Using now that H = a /a, the scale factor 
of the universe a(i), evolves with time as 



a(t) = ai (e" 1 */ 2 - l) 



2/7 



where a\ is the constant of integration. From these equa- 
tions, we can easily write the corresponding Hubble flow 
as a function of the scale factor: 



m = * 

7 



i + i- 

Ol 



-7/2' 



(53) 



Evaluating eq. (|53p at the present time (a = 1), we obtain 

7#o 



"1 



l + al' 2 



(54) 



Now utilizing equations l|50|) and (|54[) and taking into ac- 
count that the current value of the vacuum energy density 
is Aq = 3-ffgfiA (^a = 1 — Qm), we arrive at 



a 1 



3Q, 



7 - 30 r , 



2/7 



i?o(7 - 30 r , 



(55) 



Obviously, the above equation implies the following con- 
dition: 7 > 3f2 m (or n-i < 3J7a)- 

The normalized to unity, at the present epoch, scale 
factor of the universe becomes 



a(t) 



7 - 3f>, 



2/7 



(7-3fJ m )t/2 



2/7 



or 



where 



t(a) 



#0(7 - 3 ^m) 



-111 



H 



7 



7 



,-7/2 



(57) 



(58) 



It is interesting to point here that the current age of the 
universe [a — 1, E(l) = 1] is 



-In 



ff (7-3ft m ) V3tt 



7 



(59) 



We now study the conditions under which an inflection 
point exists in our past, implying an acceleration phase of 
the scale factor. This crucial period in the cosmic history 



(7/ 



3(7 - 2)On 

2( 7 - 30 m ) 



2/7 



tl 



#0(7 - 3fi r( 



-ln( 



7> 



(60) 

which implies that the condition for which an inflection 
point is present in our past is 7 > 2. 

Performing now our statistical analysis we attempt to 
put constrains on the free parameters. In particular, we 
sample 7 £ [2, 5] in steps of 0.01. Thus, the overall likeli- 
hood function peaks at fl m = 0.32±°;°2, 7 = 3.44 ± 0.02 
( 52 ) with xL( fi m,7) = 432.7 (do/=351). Using the latter 
cosmological parameters the corresponding current age 
of the universe is found to be to ~ 14.3 Gyr while the in- 
flection point is located at (ai,tj) ~ (0.48, 0.42to). This 
corresponds to a redshift Zj ~ 1.08, which is substantially 
higher than in the case of the concordance model. 

Following the notations of [5l[, we now derive the 
growth factor of fluctuations in the power series model. 
In particular, we change variables from t to a new one 
following the transformation 



y = exp(nit/2) with < y < 1 



(61) 



In this context, using equations ©, ([5H)l . (fS"3")) and (j6"T]) 
we obtain 

n\y 



Q 



Q 



7 (y-l) 2 

^1 [(6 -7)2/ -7] 
2 1 (y -1) 

"1(7 ~3)y 
2 7 (y-l) 2 ' 



(62) 



(for rii see eg 151 
We can now re- write eq.© as: 



(56) 7 2y( y _i)2^" + 2 7 ( 2; _i)(5y_ 7 ) D '_2(6- 7 )( 7 -2 2 ;)^ = 

(63) 

where prime denotes derivatives with respect to y. No- 
tice, that this variable is related with the scale factor 
as 



y=l + 7 ~ 3 ^ m a 7/2 



(64) 



We find that eq. (|63|) has a decaying solution of the form 
D x (y) = (y - l)^s)h ^ a (7-6)/2 for 7 < g. Thus, the 
corresponding growing mode of ea. (|63[) is 



D(a) = Ca (7 - 6)/2 



d:r 



where 



(2 7 -4)/ 7 



(65) 



(66) 
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with 

n rn = 1 - n A = ^ . (67) 

7 

It is interesting to mention that, the normalized Hubble 
flow [see ea. (l58l) and appendix A] can be cast as: 

E(a) = f2 A + & m a- l/2 . (68) 

It becomes clear that for 7 — > 3 we have ft m ~ Cl m . Using 
our best fit values we find Ct m ~ 0.28. Finally, in the 
limit of n\ — > (or rt2 — > 37; 7 is a constant) the current 
cosmological pattern tends to the model (see section 
4.3) as it should [A(a) ~ 37.fr 2 , E(a) ~ or 3 ^-^/ 2 and 

4.5. The A cx H model 

As we have previously mentioned, different authors 
[53l H3, [H, HI] have tried to provide some fundamental 
reason for the possibility that the vacuum term could be 
proportional to the Hubble parameter, A(o) oc H. Note, 
that this kind of cosmological model (hereafter Aps 2 ) is 
a particular case of the power series model studied in the 
previous section by setting ri2 = in eq. (|50[) : 

A = mH. (69) 

In this case, the normalized Hubble flow simply reads 
as in (|6"8"|) , but using the original (untilded) cosmological 
parameters and choosing 7 = 3: 

E(a) = Q A + n m a~ 3/2 . (70) 

One of the merits of the Aps 2 vacuum model is that it 
contains only one free parameter, n\. Obviously, this 
parameter is determined to be 

m=3H Sl A . (71) 

Therefore, since we use the prior Hq = 70.5 km/sec/Mpc 
0, the free parameter is actually Cl A (or, equivalently, 
Q m = 1 — CI a)- This single parameter is the same as 
in the standard flat ACDM model, except that the Hub- 
ble rate in ([TO]), after squaring it, compares very differ- 
ently with eq. (jTHJ) . Needless to say, this could make a 
dramatic difference when we try to fit the the combined 
SNIa+BAO+CMB data with the A PS2 model. 

Indeed, if we now marginalize the results of the pre- 
vious section over 7 = 3, the joint likelihood analysis 
provides a best fit value of f2 m = 0.34 ± 0.01, which is 
significantly larger than in the ACDM case, but with a 
poor quality fit: Xtot 

(Q ro ) ~ 513.6 for dof= 352 [p. 
This simply means that the functional form ([70)) is un- 
able to fit the observational data simultaneously at low 
and high redshifts. We confirm this point by using the 
CMB shift parameter only, finding that the correspond- 
ing likelihood function peaks at £l m ~ 0.80. This value is 



^2.5 times larger than that provided by the SNIa+BAO 
solution fl m ~ 0.32. This fact alone suggests that in spite 
of the various adduced motivations in the literature for 
the class of models A cx H, they are unfortunately unable 
to provide a quality fit of the basic cosmological data in 
all the relevant redshift ranges. 

We note that although our combined 
SNIa+BAO+CMB likelihood analysis provides a 
similar fl m value to that found in [491 ] . it has a quite 
large reduced \ 2 (— 1-46), which is in contrast with 
Carneiro et al. who concluded that the A cx H model 
was compatible with the data they used. However, 
this conclusion immediately paled after performing the 
study of the matter power spectrum of such model, 
which turned out to be highly unfavorable owing to a 
significant late-time depletion of power as compared to 
the standard ACDM model (50|. We find that this is 
borne out by the present study; in fact, the growth rate 
of galaxy clustering becomes frozen at late times for 
that model (as we will see in section [5]). The analysis 
of the fluctuation field D(a) is performed exactly as in 
section 4.4 in the limit 7^3. 

4.6. A power law A(£) model 

In this phenomenological scenario, we generalize the 
ideology of (28l . |47| in which the vacuum energy density 
evolves as A(a) cx a~ 2 . In the present work, the vacuum 
energy is taken to evolve with an arbitrary power of the 
scale factor: 

A(o) = 3 7 (3 - n)a~ n (72) 

(hereafter A„ model; see also [63j]). The corresponding 
Hubble flow as a function of the scale factor follows from 
solving eq. ^ with A given as in eg. ([72")) . Trading the 
cosmic time variable for the scale factor, the differential 
equation ^ can be seen as 

jtt2 o 

^- + ±H 2 ^3 1 (3~n)a- n - 1 . (73) 
da a 

The corresponding solution satisfying the boundary con- 
dition H(a = 1) = Ho reads as follows: 

H 2 (a) = (Hq - 3j)a~ 3 + 37a"" . (74) 

Using now the following parametrization, 

fU=l-ft m = |X , (75) 

we simply derive 

E 2 (a) = £l m a~ 3 + Q A a- n . (76) 

Evidently, if we parameterize the constant n according 
to n = 3(1 + w) then the current A„ cosmological model 
can be viewed as a classical quintessence model (Pq = 
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wpQ) with w < —1/3 (or n < 2), as far as the global 
dynamics is concerned [63j), despite the fact that the 
two models have a different equation of state parameter. 
On the other hand, utilizing eg. (fT2"|) at the present epoch 
[Ao = 37(3 — n)] and taking into account that the current 
value of the vacuum energy density is Ao = 3HqQ,a we 
obtain 



n A H 2 



and 



7 = 



SIa = 1 — tt m — 



3il A 



(77) 



(78) 



Obviously, for n — > we get fl A ~ (or f2 m ~ f2 m ) as 
we should. 

It is worthwhile to compute the evolution of the matter 
density as a function of the scale factor. From |T]), l[72|) 
and ||71J| we find 



pm{o) = 3 (ff Q — 37) a 3 + 3717a 



(79) 



We remark that, in spite of the simple form of the Hubble 
rate (|76p in terms of the formal parameters {yi m ,Q. A ) 
(see Appendix A) , the matter density is a mixture of the 
canonical a~ 3 component and the new a~™ component 
associated to the vacuum energy. This was expected from 
the fact that matter and vacuum energy are in interaction 
in this model. The physical parameters {p. m ,Q, A ) are 
related in the standard manner. Indeed, since p m o = 
Pm(a = 1), we have from (|79[) 



PmO 
3H$ 



1 



7(3-ra) 



HI 



i-n A , 



(80) 



which is consistent with the result found above for Q\. 

Comparing the A n model with the observational data 
(we sample n £ [—0.2, 2] in steps of 0.01) we find that the 



best fit values are VL r 



0.29 



+0.01 

-0.02 



(or £l m ~ 0.28) and 



n = -0.06 ±0.04 with xL(^m,n) ^ 431.2 (do/=351). 
The standard ACDM case is obtained from this model 
for n — 0, and the best fit value is indeed close to it, 
but the fact that n approaches from below (implying 
w ~ —1.02) means that the preferred situation is slightly 
tilted into the phantom domain. The inflection point 
for this model can be easily computed without resorting 
to the time dependence as follows. Requiring that the 
deceleration parameter 



aa 
• 2 



1 



dH 2 



2H 2 {a) da 



(81) 



vanishes, we find that adH 2 /da + 2H 2 = 0, and hence 
a j is the root of this equation. Using (f?4")) and (fTT)) . we 
are immediately led to 



3(2 - n)fi/ 



l/(3-n) 



(82) 



which implies that n < 2 for f2 m S (0, 1). Using our best 
fit values, we find aj ~ 0.60 (zj ~ 0.66), while the age of 
the universe is to ~ 13.6Gyr. We remark that, contrary 
to the previous models, the power law model predicts a 
transition point from deceleration to acceleration located 
at a time more recent than the concordance model. 

Finally, if we change the variables from t to a then the 
evolution of the mass density contrast (see eqf5]) becomes: 

a 2 H 2 D" + a[3H 2 + aHH + HQ]D - G{a)D = (83) 
where 



G{a) = t Y~ 2HQ - aHQ 



(84) 



Due to the fact that the best fit value n — —0.06 is 
relatively small, we can neglect from eq.(|83p the corre- 
sponding high order terms (n 2 , n 3 , n£l m £l A , etc). Thus, 
we obtain the following accurate formula to the growth 
factor (see also (63|) 



D(a) ~ a F - 



1 w- 1 



3w 2w 



1 



5 

6w ' 



f^A a 



-3w 



a 



(85) 



where w = — 1 + n/3 and -F is the hypergeometric func- 
tion. Notice that -D(a) ~ a for J1a — * 0, as expected. 



5. LINEAR GROWTH RATE FOR THE TIME 
VARYING VACUUM MODELS 

In the upper panel of figure 1 we present the growth 
factor evolution, for the different time varying vacuum 
models, presented in this work, as a function of redshift, 
D(z) [z = a^ 1 — 1]. Notice, that the growth factors are 
normalized to unity at the present time (z = 0). We 
find that for z < 0.5 the A PSl and A PS2 growth factors 
reach a plateau, implying that the matter fluctuations 
are effectively frozen, as we already advanced in section 
4.5. Also, it is obvious that the growth factors for the 
latter vacuum models are much greater with respect to 
the other 3 models, A, A^g and A„. Therefore, it is 
expected that this difference among the above cosmolog- 
ical models will affect also the predictions related with 
the formation of the cosmic structures (see next section) , 
due to the fact that they trace differently the evolution 
of the matter fluctuation field. 

We would like to end this section with a discussion on 
the evolution of the well known indicator of clustering, 
namely the growth rate (57| : 



dlnD 
dlria 



(86) 



From the known analytical form of the growth factor 
D(a) for the current vacuum models 2 , in the bottom 



The formula -jfcF (a, b, c, x) = ?fF (o + 1, 6 + 1, c+ 1, x) for dif- 
ferentiating the hypergeometric function is used to compute f(z) 
from D(z). 
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FIG. 1: Upper Panel: The evolution of the growth factor. 
The lines correspond to Apsi (solid), Arg (dashed), A„ (dot- 
dashed) and Aps 2 (dot). Note, that open triangles correspond 
to the traditional A cosmology. Bottom Panel: Comparison of 
the observed (solid circles [64|, (see Table I) and theoretical 
evolution of the growth rate of clustering f(z). We do not 
plot the growth rate of the A cosmology in order to avoid 
confusion. 



TABLE II: The reduced \ 2 values (Xmin/df) and KS proba- 
bilities comparing the growth rate of clustering between data 
and vacuum model expectations. 



Model xLu/5 Vns xlurM* < 3) Vk.s{z < 3) 



A 


0.63 


1.0 


0.10 


0.997 


A RG 


0.64 


1.0 


0.10 


0.997 


A PSl 


9.7 


0.210 


11.6 


0.107 


Aps 2 


20 


0.036 


22.9 


0.011 


A„ 


0.63 


1.0 


0.10 


0.997 



n) obtained in section 4 for each model. We then com- 
pute the corresponding consistency between models and 
data (Xmin/^ and 'Pks) and place these results in Table 
II. Obviously, in the case of A, Arc and A„ models we 
find good consistency (Xmm/5 — 0.64 and Vks — !)• 
Therefore, it is apparent that the power series Aps t and 
Aps 2 vacuum models (and specially the latter, which is 
the model where A evolves linearly with H) fail to fit the 
data {Xmin/5 ^ 20 and Vks - 0.036 for A PS2 ). In this 
framework, close to the present epoch z < 0.5, the Aps\ 
vacuum model is unable to fit the data (Xmin/5 — 9-7 
and Vks — 0.21) because the matter fluctuation field is 
frozen and thus the growth rate of clustering effectively 
overs. However, increasing the free parameter 7 by a fac- 
tor of ~ 1.2 (7 ~ 4.1) the Aps 1 model appears to fit the 
f obs (z) data (xL„/5 ^ 1.2 and V KS - 0.99). Note, that 
in the last two rows of Table II, we list the correspond- 
ing results by excluding from the statistical analysis the 
observed growth rate of clustering at z — 3 (see Table I 



TABLE I: Data of the growth rate of clustering 'fi4\. The 
correspondence of the columns is as follows: redshift, observed 
growth rate and references. 



fob 



Refs. 



0.15 0.51 ±0.11 [65,66] 

0.35 0.70 ±0.18 [67] 

0.55 0.75 ±0.18 [68] 

1.40 0.90 ±0.24 [69] 

3.00 1.46 ±0.29 [70] 



panel of Fig. 1 we display the predicted growth rate: 

d\nD 



f(z) = -(l + z)- 



dz 



(87) 



together with the observed fobs(z) [filled symbols] using 
the recent results of the 2dF and SDSS galaxies [64| . In 
Table I, we quote the precise numerical values of the data 
points, with the corresponding error bars. We compare 
the growth rate of clustering between data and models 
via a x 2 minimization and a Kolmogorov-Smirnov (KS) 
statistical test respectively. For the model predictions 
we use the best fitted valued for the parameters (7 or 



THE FORMATION AND EVOLUTION OF 
COLLAPSED STRUCTURES 



In this section we study the cluster formation processes 
by using the usual Press-Schecther formalism [71(, which 
studies the the behavior of the matter perturbations as- 
suming a Gaussian random field background, but applied 
also within the framework of the vacuum models studied 
in this work. 

We wish here to estimate, within the different vac- 
uum models, the fractional rate of cluster formation (see 
[zElzl])- In particular, these studies introduce a method- 
ology which computes the rate at which mass joins viri- 
alized halos (such as galaxy clusters), which grow from 
small initial perturbations in the universe, with matter 
fluctuations, (5, greater than a critical value S c . 

Assuming that the density fluctuation field, smoothed 
at the scale R (corresponding to a mass scale of M — 
4npR 3 /3, with p the mean background mass density of 
the Universe), is normally distributed with zero mean, 
then the probability that the field will have a value S at 
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FIG. 2: The predicted fractional rate of cluster formation as a 
function of redshift for the current cosmological models (using 
as = 0.80). The meaning of the various lines is as in Fig. 1. 



any given point in space is: 
1 



V{5,z) 



2tt(t(R,z) 



exp 



2<j 2 (R,z) 



where the variance of the Gaussian field, a 2 (R,z), is 
given by: 



a 2 {R,z) 



1 

2^2 



k 2 P(k,z)W 2 (kR)dk 



(89) 



with P{k 1 z) the power spectrum of density fluctuations 
which evolves according to P(k,z) = P(k,0)D 2 (z), with 
D(z) the growing mode of the density fluctuations evolu- 
tion, normalized such that D(0) = 1. Finally, W(kR) is 
the top-hat smoothing kernel, given in Fourier-space by: 



W(kR) 



(kRf 



[sin(fc-R) - kRcos(kR)} 



(90) 



We can now estimate what fraction of the Universe, at 
some reference redshift z and above some mass thresh- 
old M, has collapsed to form bound structures. To this 
end we need to integrate the probability function given 
by eg. ([55]) over all regions that at some prior redshift 
had overdensities which by the reference redshift have 
increased to above the critical value, 6 c {z), which in an 
Einstein-deSitter universe is ~ 1.686 and varies slightly 
for different values of f2 m (75|. Therefore, this fraction is 
given by [74j : 



T(M, z) 



V(6, z)dS 



(91) 



and performing the above integration, parametrizing the 
rms mass fluctuation amplitude at R — 8 hT x Mpc, which 
can be expressed as a function of redshift as <j(M, z) — 
<j%{z) = D{z)a%, we obtain: 



1 — erf 



V2a 8 (z) 



(92) 



Obviously the above generic form of eq. (j92]) depends 
on the choice of the background cosmology and the 
power-spectrum normalization, which we take to be the 
WMAP5 result of cr 8 ~ 0.80 g|. 

The final step is to normalize the above probability 
to give the fraction of mass in bound structures which 
have already collapsed by the epoch z, divided by the 
corresponding fraction in structures which have collapsed 
at the present epoch (z = 0), 



F{z) = T(z)/T{Q) 



(93) 



In Fig. 2 we present in a logarithmic scale the behavior 
of normalized structure formation rate as a function of 
redshift for the present vacuum models. 

In the context of the power series vacuum energy mod- 
els (see h-psx solid line and Aps 2 dot line in Fig. 2), we 
find that prior to z ~ 0.5 the cluster formation has effec- 
tively terminated due to the fact that the matter fluctu- 
ation field, D(z), effectively freezes. Also, the large am- 
plitude of the Ap^ and Aps 2 fluctuation field (Fig. 1) 
implies that in these models galaxy clusters appear to 
form earlier (z > 4) with respect to the A, Arc and 
A n (dashed lines) vacuum models. Indeed, for the latter 
cosmological models we find that galaxy clusters formed 
typically at z ~ 2. Finally, it is worth noticing that for 
a higher value of crs(> 0.80), the corresponding cluster 
formation rate moves to higher redshifts and obviously, 
the opposite situation is true for er 8 < 0.80. 



6.1. The halo abundance and its evolution 

From the previously presented Press-Schecthcr formal- 
ism (hereafter PSc), using the fraction of the universe, 
F(M,z), that has collapsed by some redshift in halos 
above some mass M, we can estimate the number density 
of halos, n(M, z), with masses with a range (M, M+5M), 
by the following: 



n(M, z)dM 



dT{M,z) p 
dM M 



dM 



(94) 



Performing the differentiation and after some algebra we 
derive the following: 



n(M, z)dM 



p / 1 da\ 



M \adM J 
diner- 1 



f P Sc{a)dM 



M dM 



f PSc {a)dM 



(95) 



where /p Sc (ct) = v r (2/n)(5 c /a) exp(-5 2 /2a 2 ). Note that 
within this approach all the mass is locked inside halos, 
according to the normalization constraint: 



-foo 



/pSc(cr)dlnCT 1 = 1 



(96) 



Although the above (Press-Schecther) formulation was 
shown to give a good rough approximation to the ex- 
pectations provided by numerical simulations, it was 
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later found to overpredict/underpredict the number of 
low/high mass halos at the present epoch (eg. [73[ and 
references therein). There is a large number of works pro- 
viding better fitting functions of /(c), which are mostly 
based on a phenomenological approach (see [77[ and ref- 
erences therein). 

From the halo mass function we can now derive an 
observable quantity which is the redshift distribution of 
clusters, Af(z), within some determined mass range, say 
Mi < M/Mq < M 2 . This can be estimated by integrat- 
ing, in mass, the expected differential halo mass function, 
n(M, z), according to: 



dV f Ah 
Af{z) = — / n(M, z)dM 



dz 



(97) 



M 1 



where dV/dz is the comoving volume clement, which in 
a flat universe takes the form: 



dV 2 dr 

-= 47 rr (z)-(z) 



(98) 



with r(z) the comoving radial distance out to redshift z: 



r{z) 



dx 

H J E(x)' 



(99) 



We use the halo mass function of Reed et al. [77[ , inte- 
grating for cluster-size halo masses, ie., between: 10 



13.4 



< 



M/Mq < 10 16 . In Fig. 3 we show the integral halo mass 
function, n(> M), for the models discussed previously 
and at two different redshifts. We use the WMAP5 nor- 
malization of the power-spectrum, ie., <7g ~ 0.8. Note 
that the classical A, the Arc an d A„ models provide 
indistinguishable halo mass functions, with those corre- 
sponding to the Ap52 and Apsi models being way off 
(we plot results only of the former model). 




10 i4 10 is 10 i8 10 i3 10 14 10 15 10" 
M/M M/M 

FIG. 3: The halo mass function at two different redshifts. The 
different models are represented by the same line-types as in 
the previous figures (which however fall on each other and 
thus are indistinguishable), but only the Aps2 model (upper 
curve) provides a significantly different n(M, z) that the other 
models. 

In Fig. 4 we show theoretically expected cluster red- 
shift distribution for the three closely resembling mod- 
els, ie., the standard A model, the Arg and A„ models 



(left panel) and the fractional difference between the first 
(constant A) and each of the other two models (right 
panel). The expected differences are small at low red- 
shifts, but become gradually larger for z > 2, reaching 
variations of up to ~ 100% at z ~ 5. To investigate 



_ 0.5 





FIG. 4: The expected cluster redshift distribution (left panel) 
and the corresponding fractional difference (right panel) of 
the Aug (upper blue curve) and A n (lower magenta curve) 
models with respect to the standard A model. 

how realistic it would be to detect such differences, we 
provide below the expectations for two realistic cluster 
surveys covering the hole sky (although a more realistic 
case would be to consider a solid angle of the order of 
~4x 10 3 square degrees): 

(a) an X-ray survey down to a flux of: /i; m = 3 x 10~ 14 
ergs s _1 cm -2 , as that expected from the future eROSITA 
X-ray satellite, and 

(b) a Sunayev-Zeldovich (SZ) survey with a limiting flux 
density at vq = 150 GHz of /„ ,iim = 5 mJy (as expected 
from the survey of the Southern Polar Telescope). 

To realize the predictions of the first survey we use the 
relation between halo mass and bolometric X-ray lumi- 
nosity, as a function of redshift, provided in [ll], ie: 



L(M,z) = 3.087 x 10 4 



M h hE{z) 



-, 1.554 



10 15 Af n 



ergs 



(100) 

The limiting halo mass that can be observed at redshift 
z is then found by inserting in the above equation the 
limiting luminosity, given by: L = 47rd|/ii m , with dh the 
luminosity distance corresponding to the redshift z. 

The predictions of the second survey can be realized 
using again the relation between limiting flux and halo 
mass from 781: 



2.592 x 10 8 mJy 



M 



1.876 



E 2 l 3 {z) 



(101) 

where dji(z) = g(l/(1 + z) 2 is the angular diameter dis- 
tance out to redshift z. 

In Fig. 5 we present the expected redshift distribution 
and the fractional difference between the different mod- 
els, similarly to Fig. 4, but now for the realistic case of 
the X-ray survey. Similarly, in Fig. 6 we present the cor- 
responding results for the case of the SZ survey. It 
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FIG. 5: The expected cluster redshift distribution (left panel) 
and the corresponding fractional difference (right panel) of 
the Aug and A n models with respect to the standard A model 
for the case of a realistic (future) X-ray survey with a flux 
limit of 10 -14 erg s _1 cm -2 . 
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FIG. 6: The expected cluster redshift distribution (left panel) 
and the corresponding fractional difference (right panel) of 
the Arc (upper blue curve) and A n (lower magenta curve) 
models with respect to the standard A model for the case of 
a realistic (future) SZ survey with a flux limit of 5 mjy. 



is evident that the imposed flux-limit together with the 
scarcity of high-mass halos at large redshifts, induces an 
abrupt decline of the Af{z) with z, especially in the case 
of the cluster X-ray survey. As can be seen from the left 
panel of Fig. 5, this fact cancels the possibility of observ- 
ing the intrinsic halo number-count differences between 
the different models (A, A^g and A„). However, the ex- 
pected SZ cluster number-counts (Fig. 6) show that we 
maybe able to detect significant differences in the red- 
shift range 2.5<, z<, 3 (at higher redshifts only a very 
small number of SZ clusters will be detected) at a level 
of ~ 6 — 12%, which translates in number count differ- 
ences, over the whole sky, of ~ 100 clusters. 



7. CONCLUSIONS 

In this paper, we have studied analytically and numeri- 
cally the overall dynamics of several time varying vacuum 
scenarios for spatially flat FLRW geometry, beyond the 
usual constant vacuum model (the standard ACDM cos- 
mology) . We wish to spell out clearly which are the basic 



assumptions and conclusions of our analysis. 

• We use various functional forms in order to parame- 
terize the vacuum energy density. In particular, we 
consider: (a) a vacuum model based on the renor- 
malization group (RG) in quantum field theory, 
A-rg = no + rt2 H 2 ; (b) a vacuum model in which A 
evolves proportional to the total energy density of 
the universe, hence A^ cx H 2 ; (c) a power series 
expansion in H up to the second order without con- 
stant term, Aps 1 = n\H + 112H 2 ; (d) a vacuum en- 
ergy density which decays proportional to the Hub- 
ble expansion rate as Aps 2 = n\H\ and (e) a vac- 
uum energy density which evolves with an arbitrary 
power of the scale factor as A„ = 37(3 — n)a~ n . 
These models have different theoretical motiva- 
tions, but not all of them are able to withstand 
the stringent experimental tests provided by the 
various observational sources. In this framework, 
we first perform a joint likelihood analysis in order 
to put tight constraints on the main cosmological 
parameters by using the current observational data 
(SNIa, BAOs and CMB shift parameter, together 
with the growth rate of galaxy clustering). Also, 
we find that the above models can accommodate a 
late time accelerated expansion. 

• In the case of the power series models Aps 1 and 
Aps 2 i we frnd that the amplitude and the shape 
of the linear density contrast are for both signifi- 
cantly different with respect to those of the A, Arc 
and A„ models. We also find that for z < 0.5 the 
matter fluctuation field of the power series mod- 
els practically freezes out, which indicates that the 
corresponding growth rate of clustering tends to 
zero. The latter result is ruled out by the observed 
growth rate of clustering measured from the optical 
galaxies. In the particular case of Aps 2 , even the fit 
of the cosmological parameters turns out to be of 
poor quality because it is unable to adjust simul- 
taneously the observational data at low and high 
redshift. In contrast, the A, Ana and A„ models 
match well the observed growth rate and they pro- 
vide good quality fits of the cosmological param- 
eters at all redshifts. We remark that the power 
law model A„ behaves effectively as a scalar field 
model with an equation of state slightly tilted in the 
phantom DE regime (w < —1). The idea that the 
models with a dynamical cosmological term A may 
behave effectively as quintessence and/or phantom 
DE has been described in general terms in the lit- 
erature [12] • 

• The particular case of the RG model for no = 
(or of the power series model for n\ =0), i.e. the 

= n 2 H 2 model, is ruled out at a high signifi- 
cance level because the best parameter fit value for 
712 is incompatible with the necessity of growing 
structure formation in this model. 
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• In the case of A, Ape and A„ vacuum scenarios, the 
large scale structures (such as galaxy clusters) form 
later (z ~ 2) with respect to those produced in the 
framework of the Aps ± and Aps 2 models (z > 4). 
Therefore, in view of the observational data, the 
former are much more favored as compared to the 
latter. 

• The expected redshift distribution of cluster-size 
halos in realistic future X-ray or SZ cluster surveys 
indicates that we will not be able to distinguish the 
closely resembling models (constant vacuum, quan- 
tum field and power-law vacuum) based on the evo- 
lution of cluster abundances using the first type of 
survey, but there is some limited hope with the sec- 
ond. 
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APPENDIX A: ALTERNATIVE VIEW OF THE 
BASIC COSMOLOGICAL EQUATIONS 

The goal in this appendix is to give the reader the op- 
portunity to appreciate the relative similarities and dif- 
ferences among the main time varying vacuum models 
(Ape Aps 1 and A„). In particular, we re-write the ba- 
sic cosmological equations in terms of a new set of formal 
cosmological parameters (f2 m , J^a)- In this new space, 
the basic structure of the ACDM Hubble flow §f§§ is 
more closely preserved, in the sense that it appears as 
the sum of fl A times a function (which can just be one 
in some case) and another term involving VL m times an- 
other function. The two functions are normalized to one 
at present, so that the sum rule f2 m + fl A = 1 is ful- 
filled. The change of basis from the physical parameters 
(O m , Q\) to the formal (O m , Q A ) ones involves additional 
parameters 7i of the current model. The two parameter 
spaces coincide only if we can set the values of the 7* such 
that the vacuum becomes static. Consider the following 
examples. 

• Ape model: If we use eq. (|37[). then the basic cosmo- 
logical equations (see section 4.2) take the following 




E 2 (a) = n A + f) m a- 3(1 - 7) . 
The scale factor at the inflection point is 

1/3(1-7) 



a 1 



(1 - 3 7 )q, 
20 A 



(A3) 



(A4) 



Obviously, in the ((l m ,(ljC) basis, the above equa- 
tions generalize more tightly those of the concor- 
dance A cosmology (see section 4.1) and reduce ex- 
actly to them for 7 = 0. 

• Aps 1 model: If we use eq. flBTl) . then the correspond- 
ing basic cosmological equations (see section 4.4) 
become: 



H{t) = n A H 



and 



a(t) 




2/ T 



e 7fi A -Ho*/2 
el Q, A H t/2 _ 1 



= 7O A ff t/2 



2/7 



or 



E(a) = A + O m a _7/2 



(A5) 



(A6) 



(A7) 



Notice that in this particular case it is E(a), rather 
than E 2 (a) , which appears decomposed as a sum of 
two terms in the new parameter space. A perfect 
analogy with the ACDM Hubble flow is not always 
possible. Indeed, in this model the vacuum can 
never coincide with that of the standard model, ex- 
cept for the trivial (and excluded) situation where 
its energy is zero. The scale factor at the inflection 
point is given by 



a 1 



(7 ~ 2R 
2^l^ 



2/7 



(A8) 



A„ model: In this case we utilize eq. (|78|) . There- 
fore, the normalized Hubble flow obeys 



E 2 (a) = n m a~ 3 + n A a~ n , 



(A9) 
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while the corresponding inflection point is 

1/(3-™) 



ai 



(2 - n)n A 



(A10) 



It is worth noting, that the current A„ cosmologi- 
cal model can be viewed as a classical quintessence 
model (Pq — wpQ, with w = — 1 + n/3), as far as 
the global dynamics is concerned. Since, however, 
n < (i.e. w < —1) is preferred by the data, in 
practice it behaves effectively as phantom DE. The 
standard ACDM cosmology is recovered from this 
model in the limit n — ► and 7 = £IaH 2 /3. 



APPENDIX C: UNIFICATION OF THE VACUUM 
MODELS 

In this appendix we examine a more general class of 
vacuum models: 



A = n + niH + n 2 H 2 , (n > 0) 



(CI) 



The time evolution equation for the Hubble flow is ob- 
tained by eq.((4|) as: 



H 



dy 



+00 -Ay 2 + my + n Q 2 



(C2) 



APPENDIX B: THE A RG GROWTH FACTOR 



where A = 3 — n 2 . If we assume that A > (or n 2 < 3) 
then corresponding general solution of eq. (|C2p is 



With the aid of the differential equation theory we 
present the growing model solution that is relevant to 
cq. (l40"|) for the RG model. Since the variable introduced 
in (f3"5| satisfies y > 1, it is possible to find the solution 
D(y) of the differential equation (HD)) in terms of the as- 
sociated Legendre functions of the second kind, Q^(y). 
The appropriate transformation reads as follows, 

n 5-3/3 I 

D(y) = (y 2 - l)-*3-Q£(2/), v=—, /x = ^-l(Bl) 



3/3 



where 



QZiv) = Cy 



-v-ii-Xi.p. 



(y'-l)*F(a,b,c,^) 



(B2) 



with 



1 + ^ + ^, b= i + ^ + C=V 



2 > 2 1 2 1 2 

C = e l ^2 



I and 



, . 1 / -.i> + j/+ 
r(«/ + f) 



(B3) 



Inserting eq. (IB2|1 into eq. (|B1[) and after some algebra we 
have 



D(y) = c^-D^y-^F[- + -, w - + -,- 2 

(B4) 



1 



1 1 1 



3 1 



D{y) = C(y 2 -l)^yF 



1 



+ 



1 3 1 



1 



3/3 2' 2' 3/3 2' - 1 

(B5) 

where we have used the well known linear transformation 
formula: 



(B6) 



F(a, b, c, x) = (1 — x) °f a, c — b, c 



x- 1 



Note, that in our formulation the parameters become: 
a= 37J + 2'k~37i' c= 37i + 2 an d x = p-- 



£T(i) 



p2 e A(p 2 - Pl )t/2 _ pi 
e A(p 2 -pi)t/2 _ I 



and 



where 



a(t) — a\ (i 



A(p 2 -pi)t/2 



2/A 



P2,l 



2A 



P2 > Pi 



(C3) 



(C4) 



(C5) 



and D = n\ + AXuq > is the discriminant. 

Obviously, the Apg 1 2 , Arc an d A models are partic- 
ular solutions of the general vacuum model. Indeed, we 
have: 

• Case 1: If no = then pi = and p 2 — ni/A. 
In this case, the basic cosmological equations [see 
eq. (|C3l) and eq. (|C4[) ] reduce to those found by ei- 
ther the Aps ± (for A = 3 — n2 > 0) or the Ap$ 2 (for 
A = 3) model respectively (see eq. lFTI) . 

• Case 2: If ri\ = then p 2 = —pi- We select the 
unknown constants such as no = SSIaHq (1 — 7) and 
n2 = 37. Therefore, the general Hubble expansion 
eq. (|C3[) , reduces to that derived by the Aug model 
(see eq|2"4")). 

• Case 3: If (111,712) = (0,0) then p 2 = —p% = 
v / 3no/3. Thus, for no = 3fl\H 2 , the Hubble ex- 



£ 0i 

pansion eq. (|C3[) . reduces to that derived by the con- 
cordance A cosmology (see eqfT6|). 



The analysis of structure formation for the model (IClj) in 
the most general case when all the coefficients rii are non- 
vanishing cannot be performed analytically. We shall 
report on this case elsewhere . 
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